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^ ' Memory and correlation effects in the interband absorption from quantum wells due to exciton-phonon 

dynamics are investigated. They are traced back to the frequency dependence and matrix character of the 

self energy arising in a 2nd Bom theory. It is found that interpeak absorption increases with respect to the 

, case in which this memory and correlation effects are neglected. 

C^ . 

^- 
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Jlj ' The interaction between excitonic and phonon modes in semiconductor nanostructures leads to both 

O I relaxation of the excitonic population as well as to pure dephasing of the optical polarization. In small 

quantum dots, where the interstate relaxation is tiny, clear signatures of pure dephasing effects have been 
measured, see e.g. the extensive work done in 1 1 1 and |2|. In this case, an accurate theoretical modeling is 
possible in the framework of an extended version of the Independent Boson Model 1 3 1. In quantum wells 
(QW), disorder-localized exciton states have significant spatial overlap and can be spectrally rather close, 
thus a combination of population relaxation and pure dephasing is expected. However, as a consequence 
of the so-called Markov approximation, the absorption spectrum of a QW results to be a superposition of 
TJ" ■ Lorentzian peaks, thus accounting for population relaxation only. We improve on this situation by going 

^^ [ beyond the Markov approximation in the solution of a 2nd Born density-matrix theory for excitons and 

acoustic phonons in disordered QW. 

We start here with the same model which was used in |4|: center-of-mass (COM) excitons are created 
j^ ■ in disorder eigenstates \a) of the QW by B^^ and are coupled to longitudinal acoustic phonons (a^ creates 

a bulk like phonon state |q)) via deformation potential interaction: 



^: n = Y,e^BiB^ + Y,hw^ a\a^ + Y.^lp{a\ + a-^)BiBp . (1) 

P^ ' a q a/3q 

' The eigenenergies Eq. are obtained from a microscopic simulation of the COM exciton in a QW with 

interfacial disorder fsl and the matrix elements t^„ ensure spatial overlap of the exciton states \a) and \P) 
and the phonon wavefunction |q}. Final output of our theory is the linear absorption spectrum a(w), which 
^ ' is connected to the Fourier transform of the exciton polarization Pa{t) = {B\ (t)) via 

cd ' 

a(a;) = ImN^ 771^^0(0') (2) 

a 

where nia are matrix elements of the dipole coupling to the light field. Both TOq, and t^^ are given expUcitly 
in 0. When deriving the equation of motion for Pa{t), a hierarchy of equations is obtained in which the 
polarization couples to expectation values of mixed exciton-phonon operators. Combining two phonon 
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operators into the phonon occupation function n{u!q) = (aJjOq) allows one to end up with a differential 
equation for Pa (t) containing a time integration over the past and cross correlations to all other exciton 

states (3 j^ a: 

{-ifidt - ea)Pc.(t) ^ m^m + n-^Y. jdt'^ap{t-t')Pp{t'). (3) 

Here we have considered an impulsive optical excitation of unit strength at t = 0. The kernel E^/j (t) 
contains 2nd powers of the exciton-phonon coupling t^^: thus, this truncation scheme is called 2nd Born 
approximation. When memory and correlation are neglected, the kernel reduces to 



'S'apit) — > 6ap5{t) 



M 



A„ + i 



M 



^^^ = — ^|iS^nn(wq)^(e,^-e„-;it^q) + (l + n(wq))(5(e„-e„ + ?iWq)] (4) 

where Aq is a polaron shift (in the following neglected because tiny) and 7^^ is the Fermi's golden rule 
scattering rate. The label M reminds that this further simplification is also called Markov approximation. 
Clearly, this leads to an absorption spectrum which is a superposition of Lorentzian peaks. However, finer 
structures are found in the spectrum a{uj) if the Markov approximation Eq. is not carried out. For 
computing them, we find convenient to Fourier transform Eq. Q, leading to the frequency-domain version 
of the kernel 



s]„^(w) = Yl h^ 



n{E)Jlp{E) 
Tiuj - iQ+ - er^ + E 



(5) 



The kernel or self energy Eq,^(ci;) depends on a coupling function J^piE) which decides which exciton 
states are coupled to the phonon bath, 

q 

In particular, J'^AE) is proportional to the overlap between |a) and |t^) and between \ri) and \(3). The 
self energy approach leads not only to Lorentzian peaks (or Zero Phonon Lines, ZPLs), but also to satellite 
absorption bands (broad bands) in a{uj). However, since 0{N^) components of J have to be computed, 
this level of approximation is computationally expensive when the number N of exciton states gets large. 
Further, it would be interesting to have a level of approximation in which the cross-correlations among 
exciton states predicted from Eq. (|3jl are neglected, but the memory effects are still accounted for. Thus, 
we consider also a diagonal (D) non-Markov approximation, which has got a numerical cost 0{N'^) but 
still reproduces both Lorentzian ZPLs and broad bands in the absorption with acceptable accuracy. In 
this case, the self energy is diagonal, i.e. Y^^Jui) — Sap'^aai'^)- The differences between the spectrum 
ai?(tL') computed using the full self energy Eq. (|5Jl and the spectrum aoi'^) obtained within the diagonal 
approximation are due to off diagonal elements of the Green's function 

[G[uj)\^p = {huj - iO+ - ea)Sal3 - ^apiuj) , (7) 

which vanish within the diagonal non-Markov approximation. Eq. is derived from Fourier transforma- 
tion of Eq. (|3} and by noticing that Pa{(^) = J2p Gai3{i-o)rap. 

For evaluating the presented theory, we consider a system of 10 exciton COM states obtained from 
diagonalization of a two-dimensional Schrodinger equation with correlated disorder potential, along the 
lines described in ||5l. The 10 eigenfunctions \a) are used for evaluating the coupling functions J'lAE) 
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Fig. 1 (Color online) (a) Interband absorption spectrum a(aj) of 10 exciton states of a 5nm GaAs/Alo.gGao.rAs 
QW at a lattice temperature T =50 K. The full line refers to the full (F) self energy computation, dash-dotted is the 
diagonal non-Markov (D) approximation, while dashes give the Markov (M) spectrum (superposition of Lorentzians). 
The 9th peak is not resolved because its oscillator strength mg is vanishing small, (b) Localization lengths A^ of the 
10 states calculated as inverse participation ratios, (c) Relative accuracy p — {ud — aF)/aF- The vertical cuts refer 
to the frequencies lja and lub used in Fig.|2| 



and, using also the eigenenergies ea, the self energy and the absorption spectra can be calculated. In Fig.[^ 
we compare the Markovian absorption (superposition of Lorentzian peaks) to the spectra obtained using 
non-Markovian self energies (both D and F version). The most striking feature is the enhancement of the 
inter-peak absorption with respect to the Markovian case. This effect is more important for states which 
are spatially less extended, as indicated by the localization lengths Aq displayed in Fig.H?. The formation 
of the broad bands is mostly due to the frequency dependence in the self energy denominator (cp. Eq. (|5})- 
This frequency dependence is a many -body effect indicating that memory in the exciton-phonon kinetics is 
not negligible. Since this phenomenon does not need to be related to relaxation among different states (set 
all indexes equal to a in Eq. (|5}), it is also said to be an effect of "pure dephasing". Minor effects present 
in the non-Markovian spectra (both D and F) are tiny renormalization of the peak positions (polaron shifts) 
and small dispersive features in the ZPLs: Both features are due to the real part of the self energy. Any 
residual difference between the D and F spectra is due to the cross-correlations predicted from Eq. (|3}- 
In Fig.[n; we show that they are responsible for deviations up to about 20% in the numerical values of 
the broad band absorption, while the ZPLs are mostly correctly accounted for The physical origin of 
the correlations is that the self energy approach accounts for scattering events \a) -^ \ri) -^ |/3), where 
a third state \ri) can "bridge" between initial and final state. Thus, it is not necessary that \a) and |/3) 
do overlap for scattering into each other Due to the t matrix elements appearing in Eq. (|6}, just \a) and 
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Fig. 2 Green's functions |Gq/3(lj)[ at frequency uj — loa (left matrix) and uj — lob (right matrix) defined by the 
cut positions in Fig.Q Logarithmic gray scale over 5 orders of magnitude; values normalized to the maximum of 
\Gai3{u}B)\- For a strictly diagonal self energy, all off diagonal elements should be zero (white). 

1 7]), and \r]) and |/3) independently, have to overlap. Conversely, within both diagonal non-Markov and 
Markov approximation only one step transitions |a) -^ |/3) are allowed in which |a) and \f3) directly 
overlap. A powerful tool for checking which states are correlated is the evaluation of the Green's function 
matrix Gafiioj). In Fig.|2]we display it for a frequency uja where D and F spectra of Fig.[^ agree, and 
another frequency ujb where large differences occur, and correlations must be important. Indeed a strong 
coupling in the subspace of states |6}, |8}, |9}, |10} is observed. The corresponding wavefunctions exhibit 
large spatial overlap and |6) has the character of "bridge-state" (not shown). 

In conclusion, we have shown that memory and correlation effects, predicted from a non-Markovian 
theory of the exciton-phonon kinetics, affect the Q W absorption spectrum. Enhanced inter-peak absorption 
is found using realistic wavefunctions for a disordered QW, which is due to memory. Furthermore, two- 
step scattering events lead to correlations in the optical response. It should be possible to demonstrate these 
predictions by means of a QW absorption measurement, which became available recently ||6l. 
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